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ABSTRACT 

Wc have made a serendipitous discovery of an enigmatic; outflow source, IRAS 05506+2414 (hereafter 
IRAS 05506), as paxt of a multi- wavelength survey of pre-planetary nebulae (PPNs). The HST optical 
and near-infrared images show a bright compact central source with a jet-like extension, and a fan-like 
spray of high- velocity (with radial velocities upto 350 km s"-'^) elongated knots which appear to emanate 
from it. These structures are possibly analogous to the near-lR "bullets" seen in the Orion nebula. 
Interferometric observations at 2.6 mm show the presence of a high- velocity CO outflow and a continuum 
source also with a faint extension, both of which are aligned with the optical jet structure. IRAS 05506 
is most likely not a PPN. We find extended NH3 (1,1) emission towards IRAS 05506; these data together 
with the combined presence of far-IR emission, H2O and OH masers, and CO and CS J=2 1 emission, 
strongly argue for a dense, dusty star-forming core associated with IRAS 05506. IRAS 05506 is probably 
an intermediate-mass or massive protostar, and the very short time-scale (200 yr) of its outflows indicates 
that it is very young. If IRAS 05506 is a massive star, then the lack of radio continuum and the late G - 
early K spectral type we find from our optical spectra implies that in this object we are witnessing the 
earliest stages of its life, while its temperature is still too low to provide sufficient UV flux for ionisation. 

Subject headings: ISM: clouds, ISM: individual objects: IRAS05506-I-2414, stars: formation, stars: 
pre-main sequence, stars: mass loss, radio lines: ISM 



1. INTRODUCTION 

We report the serendipitous discovery of an energetic 
outflow source from the object IRAS 05506-1-2414 (here- 
after IRAS 05506). IRAS 05506 was observed as part of our 
multi-wavelength survey of pre-planetary nebulae (PPNs) 
(Sahai et al. 2007), which are transition objects between 
the AGB and planetary nebula evolutionary phases in the 
lives of low and intermediate mass (1-8Mq) stars. These 
objects were selected from OH maser catalogs using the 
IRAS color criterion, ^25/ -^12 > 1-4. Subsequent and 
continuing observation of the resulting catalog of PPN can- 
didates has shown that the catalog contains PPNs almost 
exclusively, but the object we discuss here is a notable 
exception. 

In this paper, we report our HST imaging of IRAS 05506 
in the optical and near-infrared, which has revealed a 
morphology unlike that of any known PPN. We also 
present additional observations of optical, millimeter, and 
centimeter-wave lines as well as the radio continuum to- 
wards IRAS 05506, and discuss the possible nature of this 
source. The plan of our paper is as follows: in §2 we 
describe our observations, in § 3 we present the results of 
these observations, in § 4 we derive the physical properties 
of IRAS 05506 which show that it is most likely not a PPN 
but a very young protostar, and in § 5 we summarize our 
main conclusions. 



2. OBSERVATIONS 



2.1. Optical and Near-Infrared Imaging with HST 

IRAS 05506 (also AFGL5171) was imaged on UT date 
2002 October 20 (GO program 9463), by the High Reso- 
lution Camera (HRC) of the Advanced Camera for Sur- 
veys (ACS), which has a plate scale of O'.'025/pixel, using 
the F606W fllter (A = 0.59/im, AA = 0.157Aim; expo- 
sures were 2x400 s) and the F814W filter (A = O.Sl^m, 
AA = 0.166/xm; exposures were 2x200 s). All images 
were obtained with a 2-point dither, and the standard 
STScI/HST pipeline calibration has been applied to all 
data. We also obtained near-infrared images on UT date 
2004 March 23 (GO program 9801) with Camera 2 of 
the Near Infrared Camera and Multi-Object Spectrom- 
eter (NICMOS), which has a plate scale of O'.'076/pixel 
(GO program 9801), using the FllOW (A = l.l^m, 
AA = 0.38/xm; total exposure time was 544 s), F160W 
(A = 1.6/Ltm, AA = 0.28/um; total exposure time was 56s), 
and F205W (A = 2.1/Ltm, AA = 0.43/zm; total exposure 
time was 18 s) filters. For ease of comparison, the NIC- 
MOS images have been rotated to the same orientation as 
the HRC/ACS images. Numerous field stars present in the 
images have been used to achieve satisfactory registration 
between the optical and near-infrared images. 

2.2. Optical Spectroscopy 

We obtained long-sHt spectra of IRAS 05506 on 2004 
Nov 8 with the Echellete Spectrograph and Imager (ESI; 
Sheinis et al., 2002) on the 10m W.M. Keck II telescope 
at Mauna Kea (Hawaii, USA). The detector was a MIT- 
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LL CCD with 2048x4096 squared pixels of 15^m. Total 
wavelength coverage is -3900-10900A. The reciprocal dis- 
persion and the pixel angular scale range from 0.16 to 0.30 
A/pixel and from 0"120 to 0'.'168, respectively, for the ten 
echellette orders (15 to 6) of ESI, which was used in its 
echelle mode. The velocity dispersion has a nearly con- 
stant value of 11.5km s~^pixel~^ in all orders. We used 
a 0'.'5x20" slit with its long side aligned along a vector 
which passes through two of the most prominent emission 
knots seen towards this source, with an exposure time of 
600 s. Further technical details for these observations are 
described in Sanchez Contreras et al. (2008, in prep) and 
are similar to those described by Sahai et al. (2005). 

2.3. Millimeter- Wave Interferometry 

The OVRO millimeter-wave arraj|^ in its low (L) config- 
uration (for which baselines between antennas range from 
15 m to 115 m), was used to obtain snapshot observations 
of the CO J=l-0 line and the 2.6 mm continuum towards 
IRAS 05506 on 2003 Mar 17 and 24. The units of the 
digital spectral line correlator were arranged to provide 
a total bandwidth of 90 MHz (~ 234 km s~^) with a chan- 
nel spacing of 1 MHz (corresponding to ~2.6 km s^^). The 
3 mm continuum emission was observed simultaneously us- 
ing the dual-channel analog continuum correlator, which 
provided a total bandwidth of 4 GHz (after combining both 
IF bands). Total integration time on source and calibra- 
tors was ^ 11 hr. 

The calibration of the data was performed using the 
MMA software packag^ Data were gain calibrated in 
baseline-based mode using the quasar J0530-I-135, which 
was observed at regular time intervals of ~ 20 minutes be- 
fore and after our target. The quasars 3C273 and 3C84 
were used as passband and flux calibrators. 

Reconstruction of the maps from the visibilities was 
done using the Multichannel Image Reconstruction, Im- 
age Analysis and Display (MIRIAD) software. We Fourier 
transformed the measured visibilities with robust weight- 
ing (which is an optimized compromise between natural 
and uniform weighting) for CO and natural weighting for 
the continuum map for S/N optimization. After that, the 
images obtained were cleaned and maps restored. The (fit- 
ted Gaussian) clean beam for our CO (continuum) map 
has FWHM=3.9" x 2.9" (4.6" x 3.2") and is oriented 
at PA=-23.6°(-15.3°). The la noise in our CO (contin- 
uum) maps, as measured in channels with no signal, is 
0.45 (1.1) mJy beam~^. The conversion factor from CO 
(continuum) surface brightness to temperature units is 8.0 
(6.3) K per Jybeam^^. 

2.4. Centimeter Ammonia, Water and Methanol Maser 

Lines 

Observations of the H2O, NH3 and methanol lines in the 
K-band, i.e. (22-25) GHz, were obtained with the NRACF] 
Green Bank Telescope (GET) on February 28, 2006 and 
March 3, 2006. A 6x6 raster scan (step-size=13"), cen- 
tered on the J2000 position RA = 05:53:43.6, Dec = 
-1-24:14:44.5 was performed in right ascension and decli- 

^ http : //www. ovro . caltech. edu/mml 

^ MMA is written and maintained by the Caltoch's Millimeter Interfei 
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nation 4 times during these two days. The beam of the 
telescope was measured during pointing scans to be 32", 
so the raster was Nyquist sampled over a small region 
in the sky of slightly more than '--^l arcminute. A dual- 
beam receiver was used at (22-25) GHz, and total power 
measurements were made, relative to an off-source posi- 
tion of RA = 06:02:29.96, Dec = -^26: 14:44.2 (J2000). 
The GET spectrometer backend was used to record the 
total-power spectra in four different 50 MHz frequency 
bands. These four bands were set to Doppler track the 
22235.08 H2O line (616 ^23), the average of the NH3 
(1,1) and (2,2) metastable inversion lines at 23694.50 MHz 
and 23722.63 MHz, respectively, the NH3 (3,3) metastable 
inversion line at 23870.1292 MHz, and the three methanol 
lines at -24932.0 MHz (32, 1 -> 3i,2£^,42,2 -> 4i,3 
22.0 2i,i E). The 50 MHz bands provided velocity cov- 
erage of —650km s^^, centered at V^sfl^Okm s~^. Each 
50 MHz frequency band had 8192 spectral channels, pro- 
viding a spectral resolution of 6.1kHz (—0.08 km s~^). For 
each raster position, a total integration time of —16 min- 
utes was obtained, providing an rms antenna temperature 
of -25 mK. 

The VLA was used to observe the H2O line on June 14, 
2005, in the CnE configuration (beam=1.10" x 0.70" arc- 
seconds at PA=55°), with the correlator configured to pro- 
vide 127 spectral channels (Aj/= 48.8 kHz) across a band- 
width of 6.25 MHz (84 km s~^). Three velocity ranges, 
centered at Vlsr——'^^, 0; and -1-70 km s^^ (thus covering 
a total velocity range from -112 to -fll2km s-i) were 
searched by time multiplexing the observations. Each ve- 
locity setting was observed for about 10 minutes. The 
data were calibrated and imaged in the standard way us- 
ing the NRAO's Astronomical Imaging Processing System 
(AIPS) . The rms noise in the image of each spectral chan- 
nel was —15 mjy beam^^. 

2.5. Radio Continuum Observations 

We observed the radio continuum at 8.4 GHz (3.6 cm) 
in IRAS 05506 on June 11, 2005, using the VLA in the 
CnE configuration. Approximately 13 minutes were spent 
on source. The data was calibrated in the usual man- 
ner, using AIPS. A continuum image was made using the 
calibrated u, v data; the synthesized beam, using natural 
weighting, was 5"1 x I'.'l at a position angle of 56°. The 
rms noise in the image was —50 ^Jy beam~^, similar to 
the expected theoretical noise. 

We also examined the region around the position of 
IRAS 05506 in the NRAO VLA Sky Survey (NVSS; Con- 
don et al. 1998). The NVSS observed at a frequency 
of 1.4GHz (21cm), and typically had noise levels in less- 
confused regions of about 400 /xJy beam~^. This was ap- 
proximately the noise in the field containing the position 
of IRAS 05506. 

3. RESULTS 

3.1. Optical and Near-Infrared Imaging 

The HST F606W image (FigjT}i) shows a compact 
source {Sa), and a fan-like spray ofcompact nebulous fea- 
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tures Kl — K9 that are separated from Sa, but appear 
to emanate radially from it within a range of position an- 
gles (121°-155°). Examination of the most prominent of 
these knots {Kl — KA) shows that each of them is a highly 
coUimated structure. Most of these knots are seen in the 
other filters used for the HST imaging as well. A few 
more knots (fclO — kl2), not seen in the optical images, 
become visible in the FllOW and F160W images in the 
PA range 85°-121°(Fig.[2]Q at locations closer to Sa, pre- 
sumably because these are more highly extincted than the 
Kl — KA knots. Compact knots like fclO can be distin- 
guished from faint stars in the NICMOS images (labelled 
fsl — fsA in Fig.[2| because, relative to their surroundings, 
the knots are brightest in the FllOW image (i.e., at the 
shortest IR wavelength), whereas the stars are brightest 
in the F205W images (i.e., at the longest IR wavelength). 
Significant filamentary structures can be seen all around 
Sa. The most prominent of these lie in the western hemi- 
sphere around Sa, and several of these appear to be ra- 
dially directed away from the latter (dotted arcs in Fig[2] 
delineate angular wedges where these structures are most 
clearly seen). Two isolated compact nebulosities are also 
seen in this hemisphere, fclB and fcl4 - the latter has an 
inverted Y-shape, and may represent two (or more) ad- 
joining elongated knots. Knot kl3 appears to lie on the 
periphery of an extended faint structure (Tube), described 
below. 

A second compact source {Sb), comprised of a close stel- 
lar pair with a separation of 0'.'12 (inset, panel b), is seen 
5'.'64 to the west (along PA=— 97.6°) of the main source. 
Sa has a cometary structure, with a bright peak at the 
northeast end Sa-pk, and a jet-like extension Sa-jet to- 
wards the southwest, along PA~230° (Fig|3]), which is 
nearly perpendicular to the average direction of the knot 
spray seen in the F606W image. 

At the location of Sa-pk, the NICMOS FllOW (Fig{l]:), 
F160W and F205W (not shown) images show a point 
source of increasing brightness with wavelength. For this 
point source, we derive fluxes of 6.5, 54, and 190 Jy, re- 
spectively, at 1.1, 1.6 and 2.1^m from the FllOW, F160W 
and F205W images, using the photometric calibration in- 
formation provided by STScf^ The l.l^m flux is less re- 
liable than the 1.6 and 2.1/im fluxes because the star is 
relatively less bright compared to the surrounding nebu- 
losity. The lack of PSF structure towards Sa-pk in the 
optical images indicates that the central star is not being 
seen directly at optical wavelengths. Archival 2MASS im- 
ages of IRAS 05506 and vicinity (Fig.l4| show that Sa lies 
atop extended diffuse nebulosity which can be easily seen 
extending as far as ~ 20" — 25" towards the northwest. 

The position of Sa-pk as measured from the F606W 
HST image is (J2000) RA = 05:53:43.55, Dec =24:14:44.0. 
The 2MASS All-Sky Catalog of Point Sources (Cutri et 
al. 2003) shows only one source 

(05534356-^24144470 
within a ~ 4" radius of Sa, and since the uncertainty in 
the absolute HST astrometry is < 1", the latter represents 
the near-IR counterpart of Sa. The next nearest 2MASS 
source is 05534318-1-2414441, located - 5" to the west of 



Sa, and represents the near-IR counterpart of Sb. 

We note the presence of two intriguing features in the 
vicinity of Sa and Sb, in the FllOW and F160W NICMOS 
images (Fig.e. One is a tube-like, limb-brightened struc- 
ture (labelledas "Tube^' in the figure) that appears to orig- 
inate from near Sa and extends towards Sb. The second is 
a roughly linear structure, L, joining two locally bright re- 
gions {Bl and B2); the vector joining these regions passes 
through Sa. Bl is bright, relatively compact and appears 
to be located on the limb brightened periphery of the Tube 
feature, whereas B2 is much fainter, larger in size and dif- 
fuse and part of an extended nebulosity aligned roughly 
east-west. In the F205W image, since Sa is much brighter 
relative to the surrounding nebulosity, only those portions 
of the Tube and L features which are relatively more dis- 
tant from Sa can be clearly seen. The similarity of the fea- 
tures in the near-IR images indicates that these are dusty 
structures seen in scattered ambient starlight. We spec- 
ulate that the Tube might be a dense-walled cylinder of 
molecular gas and dust and may represent the wake pro- 
duced by the motion of Sb through the dense interstellar 
cloud detected in our NH3 mapping. The feature L might 
also be a wake produced by the motion of Sa, or may rep- 
resent a highly-coUimated ejection from Sa. Near-infrared 
emission-line imaging and spectroscopy is needed to better 
understand the nature of the Tube and L structures. 

3.2. Optical Spectroscopy 

The slit was aligned along a vector which passes through 
the knots Kl, K2 and Sa. The optical spectra of knots 
Kl, K2 show high- velocity Ha and forbidden line emis- 
sion such as, e.g., [NII]AA6548,6583A, [SII]AA6716, 673lA, 
[OI]AA6300,6364A, [Ca II]AA7291.5,7323.9A and many 
lines of [Fell] . The position- velocity structure and wide ve- 
locity widths (~350km s~^) seen clearly in the strongest 
of these lines (Ha, |NII], [SII]) are characteristic of bow- 
shock emission (Fig|6| resulting from the interaction of a 
knotty, high- velocity outflow (or outflows) with dense, am- 
bient material. The knots only show emission lines and no 
continuum, similar to Herbig-Haro objects. High-velocity 
emission is also seen at the location of Sa, most likely as- 
sociated with Sa-jet. In each case we detect blue-shifted 
emission only, which can be interpreted in terms of the 
near side of a bipolar outflow in which the far (i.e. red- 
shifted) side, if present, is hidden from our view due to 
local extinction. Unlike the knots, Sa also shows a very 
red optical continuum. 

An analysis of the spectrum (Sanchez Contreras et al. 
2008) reveals numerous metallic absorption lines due to Ti 
I, Fe I, and Ca II in the - 8400-8800A wavelength region 
where the continuum from the central star is less extincted 
and thus detected with much higher signal-to-noise than 
at shorter wavelengths. In particular, the Ca II triplet 
(8498, 8542, 8662A) is quite strong. The Paschen lines, in 
the same window, are noticeably absent - Sanchez Contr- 
eras et al. (2008) conclude that the most likely spectral 
type of the star is consistent with late G to early K. 
Long-sht optical spectroscopy of IRAS 05506 was re- 

^ in this figure, we have retained the original orientation of the NICMOS images (which is different from the optical images taken with the 
ACS) in ord er to avoid any lo ss of fidelity which would result from image rotation 
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ported by Manchado et al. (1990) with a 1'.' 5- wide slit 
oriented along PA=80°; they found an emission-line spec- 
trum from an east "lobe" , but only continuum from a west 
"lobe"; the lobes being separated by 5". Our HST images 
show that, if the "east" lobe is source in that study, as 
is likely, then the west "lobe" is source Sb. 

3.3. Millimeter- Wave Interferometry 

The CO line profile is centered at Visr — 7.6km s~^, 
and has a core with a width (FWHM) of ^15 km s""'^ and 
weak, broad wings which extend at > 50 km s~^ on ei- 
ther side of the line center (Fig.[7]). The core emission 
is compact (< 5"in size) (Fig.lsl). A lOmJy continuum 
source (at A = 2.6 mm) with an unresolved core and a 
faint extension along PA~230°, consistent with the PA of 
the optical jet feature, Sa-jet, is also seen at the same loca- 
tion. A plot of the emission in the extreme red- {Visr 49.4 
to 119.6km s-i) and blue- {Visr - 36.4 to -106.6km s'^) 
wings (Fig.|8| shows that the outflow is directed roughly 
along the same PA as Sa-jet and the extension in the con- 
tinuum source. This high- velocity bipolar outflow is thus 
directed along an axis roughly orthogonal to the average 
direction of high- v eloc ity ejections represented by the op- 
tical knot spray (§3.2|. 



3.4. Ammonia, Water, Methanol and Radio Continuum 
Observations 

With the VLA, we detected a single, spatially- 
unresolved H2O emission feature (across three spectral 
channels) at (J2000) RA = 05:53:43.54, Dec = 24:14:45.2, 
at an LSR velocity of 5.9 km s^^ and with a peak flux 
density of 145 mJy beam~^ (small circle in Figfla). The 
la error in the H2O maser position is < 0.1". Given the 
astrometric accuracy of the 2MASS catalog (< 0.1"), we 
conclude that there is a small, but significant separation 
between the locations of the H2O maser peak and Sa of 
0.57". This small separation, and the similarity in the H2O 
maser radial velocity and the systemic velocity strongly 
suggests that latter is associated with Sa. 

With the GBT, we detected emission in the H2O and 
the NH3 (1,1) and (2,2) lines. In H2O, two velocity fea- 
tures are detected; one at 6.16km s^^ and another at 
Vlsr = 11.15km s"^ (Fig.|9]). The FWHM widths (un- 
certainties) of these features, derived from Gaussian fits 
are 0.51 (0.04) and 0.56 (0.01) km s~^ Both velocity fea- 
tures are spatially unresolved. These, as well as the NH3 
1,1 brightness distribution, peak at the same raster posi- 
tion, which lies ~8"4 northwest of the VLA H2O maser 
position. But given the pointing uncertainties of the GBT 
and its large beamwidth, this offset is not significant. The 
H2O 6.16 km s~^ feature is most likely the same as that 
detected with the VLA (hereafter 6km feature). The 
integrated fiuxes (errors) of the H2O maser features for the 
2 epochs, June 14, 2005 and Feb 28, 2006, are as follows: 
0.163 (0.017) Jy km s-i and 1.03 (0.056) Jykms^^ for 
the 6km s~^ feature and <0.051 Jykm s~^ (3cr) and 12.25 
(0.075) Jy km s~^ for the 11km s^^ feature. Thus, over 
the 1.3 yr period between the two epochs, both features 
increased dramatically in strength (the 11km s~^ feature 
by a much larger factor), indicating that they are physi- 
cally associated and/or the variations represent a common 
response to the same stimulus (e.g. passage of a shock 



wave). 

The NH3 lines are surprisingly narrow: the FWHM 
of the higher signal-to-noise 1,1 line is 0.96 km s^^ (dis- 
cussed further in s|5]). The NH3 emission line profile in 
the main hyperfine component is centered at V^sR = 
6.21km s~^ (Fig. 10 1. The NH3 emission is extended to- 
wards the east and northwest; the region of emission is 
likely slightly larger than 1' (the raster was not large 
enough to measure the full extent of the ammonia emis- 
sion). We detect weak ammonia (2,2) emission (only in 
the main hyperfine component, and after smoothing the 
spectrum by a factor of four) , at the position of the peak 
of the (1,1) emission. No methanol lines were detected. 

Our 8 GHz continuum observations show no continuum 
emission (3ct upper limit of 150 /iJy) within a field of 
approximately 40" diameter centered on the position of 
IRAS 05506. From the 1.4 GHz NVSS all sky survey data, 
we obtain an upper limit of 1200 fi3y. 

4. PHYSICAL PROPERTIES OF IRAS 05506 

The steep increase in Sa-pk's flux from 0.6 to 
2fj,m strongly indicates that it is a highly extincted stel- 
lar source which most likely is driving the energetic out- 
flow responsible for the whole spray of knots. Consis- 
tent with this interpretation, we flnd a bright source 
in the MSX (Midcourse Space Experiment: Egan et al. 
1999) catalog located at the position of Sa-pk, i.e., at 
(J2000)RA=05:53:43.55, DEC=24:14:44.0. In this sec- 
tion, we derive the physical properties of the central source 
and its environment, many of which depend on the dis- 
tance to the source. We estimate a formal kinematic 
distance of 2.8 kpc to IRAS 05506, using the formulation 
given by Goodrich (1991). We have taken the average 
of the source radial velocities from our ammonia and CO 
data, i.e., 6.9 km s~^ as the systemic velocity, and assumed 
Galactic parameters i?o = 8 kpc, O ~ 220 km s~^ and 
circular rotation. But given the galactic longitude of 
1/7 = 185°, and allowing a typical molecular cloud veloc- 
ity dispersion of ^^lOkm s~^, the possible distance range 
is not well constrained. In this paper, we adopt a fidu- 
cial distance, D=2.8kpc, and scale all of our estimates 
with D. At this distance, the bolometric luminosity of 
IRAS 05506 derived from integrating over its (very red) 
SED, is, L ~ 51OOL0. The integration spans the full range 
of wavelengths covered by ISO archival spectra of this ob- 
ject, i.e., 2.4 — 198/im. The uncertainty in estimating L 
due to fiux beyond the wavelength limits of the SED is 
<5%. 

Given the spectral type (late G - early K, i.e. Teff ~ 
5000 — 4400 K) of the central star inferred from our opti- 
cal spectroscopy, and the above luminosity, we can locate 
IRAS 05506 on the HR diagram and compare with evolu- 
tionary models. Inspection of the pre-main-sequence stel- 
lar evolutionary models of Bernasconi & Maeder (1996) 
shows that the above values of T^f / and L lie between the 
tracks for stars with main-sequence masses of 9 and 15M0, 
and at the low temperature end of these tracks - thus close 
to the epoch where deuterium-burning is initiated. 

4.1. The Optical Knot Spray: Interstellar Bullets 

We propose that the knots are Herbig-Haro objects 
based on their forbidden-line emission; they are seen most 
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prominently in the F606W image compared to the im- 
ages at longer wavelengths because this filter covers the 
Ha and [Nil] lines. The [SII]A 6716/A 6731 ratio, R[sii], 
is about 0.6 (0.5) for the low (high)-velocity emission in 
the knots, implying electron densities of about 3300 (8000) 
cm~^ assuming an electron temperature of 10^ K (Canto 
et al. 1980). Using the average value of R[sii] for knot 
Kl (0.54), we derive Ue — 6000 cm^'^, and approximat- 
ing Kl's shape with a prolate ellipsoid having major (mi- 
nor) axis=0'.'29 (0'.'15), we find, assuming rie = uh, that 
its mass is, Mki = 2 x 10^^(0/2.8 kpc)3M0. The knots 
K2 — KA are somewhat fainter (by factors of 2-3) than 
K\, and therefore may have correspondingly lower, but 
comparable, masses. 

The morphology and kinematics of the fan-like spray 
of elongated knots in IRAS 05506 are similar to the in- 
terstellar bullets seen in H2(2.12/xm) and [FeII](1.65/im) 
emission towards the Orion molecular cloud (e.g., Allen & 
Burton 1993). Allen & Burton (1993) derive a minimum 
mass of ~ 1O~^M0 for the Orion bullets, comparable to 
our value for knot K\. The opening angle of the bullet 
spray in IRAS 05506 is about 45° (70°), measured from 
knot KS) to fclO (fcl2), hence somewhat smaller than the 
almost 90° opening angle covered by the Peak 1 region 
(located northwest of BN and Source I) shown in Nissen 
et al. (2007). The Orion bullets are typically 2" - 4" 
in angular extent (Tedds, Brand & Burton 1999); if we 
assume that the physical size of the IRAS 05506 bullets 
is similar, then their smaller angular sizes suggest that 
IRAS 05506 is about a factor 10 more distant than Orion, 
i.e., 4.5 kpc. The total (angular) radial extent covered by 
the bullets in the Peak 1 Orion region is about 35" — 40", 
i.e., a factor 6-7 times larger than that measured for the 
IRAS 05506 bullets region - consistent with a larger dis- 
tance to the latter by a similar factor, if we assume sim- 
ilar physical extents for the bullet regions of these two 
sources. The 6-10 times larger distance to IRAS 05506 
could also account naturally for the fact that the num- 
ber of high- velocity knots in IRAS 05506 is much smaller 
than that found toward Orion, since we presumably detect 
only the most extreme (i.e., brightest and largest) mem- 
bers of a larger collection of knots. Similarly, it may be 
that the opening angle of the IRAS 05506 "bullet" region 
is intrinsically comparable to that found for Orion, but our 
current observations lack the sensitivity needed to detect 
lower brightness knots at wider opening angles because of 
IRAS05506's much larger distance. Furthermore, if the 
filamentary structures seen in the FllOW image (Fig.[2|, 
together with knots /cl3 and klA represent the unresolved 
or partially resolved counterpart to the bullet spray seen to 
the east of Sa, then the similarity to the Orion bullet spray 
becomes even more pronounced. For knot A'l, which is the 
most distant knot from Sa, we estimate an expansion time- 
scale of 200 (D/2.8kpc) yr , under the assumption that the 
FWZI of the velocity-profile of a bow-shock is equal to the 
intrinsic speed of the bullet producing the bow shock (Har- 
tigan, Raymond & Hartmann 1987). For knot K2, which 
is closer to Sa than K\, the expansion speed, inferred as 
above, is smaller (^300 km s~^), and we find an expansion 
time-scale of about 170 (D/2.8kpc) yr. We don't know if 
all the knots were ejected quasi-simultaneously over some 

12 PiUai et al. (2006) find values in the range (0.7-10) X 10~* 



modest range of velocities from the central source, or re- 
sult from an episodic outflow which ejects high-speed blobs 
intermittently over a wide range of position angles. 

In addition, we have examined the literature on a 
few other examples of possibly similar systems with 
high-velocity knots. For example, the high-mass YSO 
IRAS 18151-1208 contains knots spread along two, roughly 
orthogonal, axes, and are related to two well-defined coUi- 
matcd outflows (Davis et al. 2004). The outflows IIH212 
and HH211, which arise from low-mass protostars, also 
have very highly coUimated outflows (Lee et al. 2006, 
2007). Thus these three YSOs are morphologically quite 
dissimilar to the knotty ejections from IRAS 05506, which 
have an opening angle of at least 44° . The high-mass YSO 
IRAS 05358-1-3843, displays a fan of knots seen in near-IR 
H2 emission, and for this source, Kumar et al. (2002) make 
a similar analogy to the Orion bullets. 

The optical jet, Sa-jet, is much bluer than the central 
source, Sa-pk, presumably because like the knots, it is 
dominated by optical line emission. The rather precise 
alignment of the position angle of this jet, which is seen 
out to a radius of 0'.'25 from Sa-pk, to that of the much 
more elongated (6", giving a deconvolved length of about 
4'.'5) feature seen in the millimeter continuum image, is 
remarkable. Is this outflow different in nature than that 
which produced the bright knots Kl — KA, since no such 
knots appear to be associated with it? It may be that such 
knots are present, but not visible due to local extinction. 
The average large-scale (i.e., over a_0.5' area) extinction 
derived from our NH3 data (see 



4.2 below) is A\ 



and it is likely that the extinction increases steeply towards 
Sa. Indeed, we find that the line-of-sight extinction to Sa 
is Av ^ 25 by comparing each of our observed 1.6 and 
2.2/im fluxes of the central star with the values expected 
based on its bolometric luminosity and effective tempera- 
ture. We have used the tabulated average extinction curve 
in Whittet (1992) to derive the visual extinction from our 
near-IR values. A faint patch of nebulosity is seen in the 
F606W image at an offset of about 3'.'5 along the PA of the 
optical jet, but it does not have a well delineated structure 
like the Kl — KA knots and is similar to other faint diffuse 
nebulosities seen in the vicinity of Sa, so it is not clear 
whether or not it is directly associated with the jet. 

4.2. The Gas and Dust Environment 

The physical properties of the molecular gas associated 
with, and/or in the vicinity of IRAS 05506 on an arc- 
minute angular scale, can be inferred from the ammonia 
data. Using Gaussian fits to the ammonia line profiles to 
estimate the peak intensities, and the formalism given by 
Wilson, Gaume, & Johnston (1993), Ungerechts, Walms- 
ley, & Winnewiser (1986), and Ho & Townes (1983), we 
estimate the opacity of the (1,1) transition (r '-^ 0.7 for 
the main line) and the rotational temperature between 
the (2,2) and (1,1) transitions {Trot = 12.5 ± 2.2 K). The 
rotational temperature is a lower limit to the kinetic tem- 
perature, Tfci„; for these low values of T^ot, T/jj„ ~Trot 
(Morris et al. 1973, Danby et al. 1988). Assuming a beam 
filling factor of unity, the ammonia column density toward 
the peak of the emission is ^nh^ ~ 1.1 x 10^'*cm~^. Esti- 
mating the region of ammonia emission to have a diameter 
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of 1', and adopting a fractional NH3/H2 abundance ratio 
(by numberf^ /{NHa) = 4 x 10-^ we obtain a molec- 
ular hydrogen mass in this region of 11 (D/2.8kpc)^Af0. 
In this calculation, we have accounted for the lower emis- 
sion at the outer regions of the source, by assuming that 
the derived mass is a factor 2 lower than the value that 
would result from assuming uniform emission within the 1' 
region. We estimate that the large-scale visual extinction 
towards the center of the cloud is. Ay ^ 3, using a ratio 
N{H2)/Av = 0.94 X 10^^cm~2 mag"^ (Frerking, Langer 
& Wilson 1982). 

The virial mass of the cloud (assuming an r^^ density 
distribution) is, M^ir = 188(R/pc) (AVi/2/km s-^)^ = 
73(D/2.8kpc)(0/O.5')Mo, where is the cloud radius in 
angular units, and AVi/2=0.96km s^^. For an r^^ den- 
sity cloud, we get M„„=48(D/2.8kpc)(6l/0.5') Mq. The 
factor 4 (or more) discrepancy between the virial cloud 
mass and that derived from the ammonia data may easily 
be attributed to the large uncertainties in our assumed 
value of f{NH^) and source distance, and/or to non- 
gravitational linewidths, such as unresolved outflow mo- 
tion. 

We have very recently (Nov 2007) mapped the ^^CO 
J=l-0 emission in a 36' x 18' (Aa x Ai5) region around 
IRAS 05506 using the ARO 12-m millimeter-wave tele- 
scope, and find the presence of a molecular cloud which 
extends over at least 30' East- West and 12' North-South. 
Observations of the ^'^CO and C^^O J=l-0 emission to- 
wards the central region show that the ^^CO line is op- 
tically thick, both in the line core and wings. We will 
discuss these data in detail in a future publication; our 
preliminary analysis shows that the extended cloud which 
we detect in ^^CO J=l-0 emission is easily many hundreds 
of solar masses. 

The molecular cloud detected in NH3 emission is pre- 
sumably associated with the extended diffuse nebulosity 
which is seen in the 2MASS images. The nebulosity is 
most likely due to scattered light from dust. In order to 
estimate the dust mass associated with this cloud, we have 
retrieved and re-calibrated archival ISO SWS and LWS 
spectra of IRAS 05506 covering the 3-200 /zm wavelength 
range (Fig 111. We find that it shows a very red con- 
tinuum, a deep 10 /im amorphous silicate absorption fea- 
ture, and emission in numerous far-infrared CO, OH and 
H2O lines, as well as the 01 (at 63 and 146 fj,m) and CII 
(158 Atm) lines (Sahai & Sarkar 2008, in prep). The far-IR 
line spectrum of IRAS 05506 is very similar to that seen 
towards the low-mass protostar IRAS 16293-2422 (Cecca- 
relli et al. 1998). A multi-component dust model (i.e., 
with "cool" , "warm" , and "hot" components, Sahai et al. 
1991) assuming a power-law (A^^, withp = 1.5) dust emis- 
sivity with a value of 150 cm^g"^ (per unit dust mass) 
at 60 /im (Jura 1986), provides a modest fit to the SED, 
and demonstrates that the bulk of the mass is in cool 
dust. We find that the cool component has 0.26 Mq of 
dust at 27 K, the warm component has 0.22 x 1O~^M0 of 
dust at 79 K; both masses scale as (D/2.8kpc)^. These 
imply a total mass of 26 (assuming a gas-to-dust ra- 
tio, gd = 100), which compares reasonably well with that 
derived from our ammonia observations. Note that the 
above modelling can provide only a rough estimate of the 



mass because the dust is expected to be distributed over 
a range of temperatures, depending on the distance from 
the heating source(s). Sahai et al. (2007) have compared 
the results of such modelling with those from more accu- 
rate modelling (using a dust radiative transfer code) for 
pre-planetary nebulae, and find that the simple modelling 
gives dust masses that are lower by factors of 0.6-0.8. 

For the mm-continuum source associated with Sa, we 
derive a dust mass of (0.022-0.064) (D/2.8 kpc)2M0, as- 
suming 79-27 K for the temperature of the dust, and that 
p = 1.5 is valid over the A = 60 /xm to 2.6 mm wavelength 
range. Thus the total mass associated with within a 
~4"-size region is 2.2-6.4 M© (assuming gd as above). The 
total molecular mass derived from the CO fiux detected 
with the OVRO beam, i.e., within a ~3'.'4-size region, is 
0.26 M0 assuming a fractional CO/H2 abundance ratio, 
/(CO) = 10-^ and taking T^ot=27K0 Because the CO 
is optically-thick, this value represents a lower limit. 

The inadequacies of the SED fit from our simple model 
indicate that the source has dust over a wide range of 
temperatures. Sahai & Sarkar (2007) have therefore used 
the online tool provided by Robitaille et al. (2007), which 
provides least-squares fits of pre-computed disk-envelope 
models of young stellar objects (YSOs) to user-defined 
SEDs, in order to fit IRAS05506's SED. Given an input 
distance range of 0.1-5 kpc, the fitting tool finds that the 
first three best-fit models to the SED of IRAS 05506 have, 
respectively, distances of D=2.2, 2.6, 3.0 kpc (implying 
bolometric luminosities in the range (3150-5850) ^0), cen- 
tral star masses of M=8.5, 10.6, II.2M0, and Te//=4190, 
4430, 4280 K. These Teff values are consistent with the 
spectral type of the central star inferred from our optical 
spectroscopy. 

4.3. Molecular Outflow: Mass and Dynamics 

We estimate the mass (Mtip), mass-outflow rate (Mtip), 
scalar momentum (Psc), mechanical energy {Emech) and 
the mechanical power (Lmech) in the high-velocity bipo- 
lar outflow using the integrated CO J=l-0 line profile 
from our OVRO data. The derived values for Mmp, Psc, 
and Emech scale as (D/2.8kpc)^, whereas Aff,jp and L^ech 
scale as D/2.8kpc. The scalar momentum and the me- 
chanical energy are computed using the methodology de- 
scribed in Bujarrabal et al. (2001). As before, we as- 
sume Tj.o(=27K, f{CO) — 10^^, and the emission to 
be optically thin - since the CO emission is optically 
thick, our derived values are lower limits. In this anal- 
ysis, we use the J=l-0 emission flux over the Visr velocity 
ranges [10.4 : 119.6] km s^^ for the red component and 
[2.6 : —106.6] km s^^ for the blue component of the out- 
flow shown in Fig.|8] We have conservatively excluded 
emission in the central three 1 MHz channels which cover 
the line core. We infer an expansion time-scale for the out- 
flow of 185 yr, by dividing the linear separation between 
the red and blue outflow components by the mean sep- 
aration velocity, 156 km s~^ (derived from the difference 
between the average velocities for the red and blue com- 
ponents in Fig.js]) and assuming an intermediate inclina- 
tion (45°) of the outflow axis relative to the line-of-sight. 
We find Mmp > O.ISMq, Mup > 9.7 x IQ-'^Mq yr"!, and 
Psc > 8.6 X 10^^ gcms-i (4.3M0km s^^). The mechan- 



if the excitation temperature is higher, the derived mass will be larger; e.g., for Trot=79K, the derived mass is a factor 2.5 larger 
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ical energy in the bipolar outflow is, Emech > 4.6 x 10''^ 
erg; dividing this by the expansion time-scale, wc obtain 
Lmech > 7.9 X 10^^ erg s~^. Hence the ratio of the mechan- 
ical power to the stellar luminosity is, Lmech I L > 0.04, 
which is significantly larger than the vahics typically found 
for YSO outflows (IQ-^ - 10-2; Lada 1985). The out- 
flow cannot be driven by radiation pressure: the expansion 
time-scale, ~185 yr, is much smaller than the time required 
by radiation pressure to accelerate the observed bipolar 
outflow to its current speed, trad = Pscl{L/c) 4 x 10^ 

yr- 

5. DISCUSSION 

The morphology of IRAS 05506 docs not resemble that 
of any known PPN. Although its luminosity (at a distance 
of 2.8 kpc) is commensurate with that expected for PPNs, 
the circumstellar mass is signiflcantly larger than that typ- 
ically found in PPNs (< IM©, e.g. Sahai et al. 2007). 
Zinchcnko, Pirogov & Toriscva (1998) include IRAS 05506 
(labelled by its Galactic coordinates, as G 184.96-0.85) in 
their study of dense molecular cores in regions of massive 
star formation, and find CO J=l-0, CS and C-'^'^S J=2-l 
emission towards this object using the Onsala 20-m tele- 
scope (bcamsizcs 35" and 40"). The radial velocity of the 
CS line, 6.2 km s~^ Visr, agrees with that found from our 
NH3 observations, so there can be little doubt that the 
CS emission is probing the same molecular cloud core as 
NH3. Our recent mapping of CO .1=1 emission towards 
IRAS 05506 with the ARO 12-m dish shows the presence 
of a very extended molecular cloud with a total mass of at 
least several hundred solar masses. 

IRAS 05506 is, in all likelihood, not a PPN. The 1612 
MHz OH maser line profile in IRAS 05506, observed with 
the Arecibo dish (Chengalur et al. 1993), does not 
show the typical double-peaked shape seen towards most 
AGB stars or PPNs. The main hues of OH at 1665 
and 1667 MHz were detected by Lewis (1997) at Arecibo; 
the emission covers about 15km s~^ in the 1612 MHz 
line and ~40-45km s~^ in the main lines. The com- 
bined presence of maser (H2O and OH), NH3, CO, CS 
& C^'*S and relatively strong far-IR emission, strongly ar- 
gues for a dense, star-forming core being associated with 
IRAS 05506. Palagi et al. (1993) also classify IRAS 05506 
as a star-formation region, on the basis of a principal 
component analysis of the IRAS colors of all H2O maser 
sources in the Arcetri Catalog (Comoretto et al. 1990) 
which had IRAS PSC counterparts. The very short ex- 
pansion time-scale of the optical knots and the fast molec- 
ular outflow indicates that the source Sa is a very young 
protostar. 

The presence of OH maser emission argues for 
IRAS 05506 not being a low-mass YSO, since there are 
no known OH masers toward low-mass YSOs. There have 
been surveys for OH maser emission from, e.g., T Tauri 
stars (Chandler et al. , private communication) , but no de- 
tections have been made and we are unaware of any pub- 
lished detections of OH masers toward low-mass YSOs. 
Given our calculated luminosity at D=2.8kpc, our dis- 
tance estimate would have to be lower by at least a factor 
4 for the luminosity to be in the range for low-mass YSOs. 

Is IRAS 05506 an intermediate-mass or massive YSO? 
In a search of the literature on OH masers towards star- 



forming regions, we found only one region (L1287) which 
does not appear to be forming massive stars and which 
shows an OH maser (Wouterloot, Brand, & Fiegle 1993). 
This region contains the FU Ori-type objects RNO IB, 
RNO IC and an IRAS source (IRAS 00338+6312) whose 
luminosity indicates it to be an intermediate-mass YSO. 
So there is a small probability that IRAS 05506 is associ- 
ated with an intermediate-mass star (Quanz et al. 2007 
and references therein). Further, the total molecular cloud 
mass derived from the ammonia data is rather small for a 
massive star, because the implied efficiency of star forma- 
tion would have to be signiflcantly higher than is usually 
assumed to be the case. Note, however, that this issue is 
a lesser concern if the source were at a distance exceed- 
ing 5 kpc, because of the scaling of the cloud mass with 
the square of the distance. Also, the extended molecular 
cloud towards IRAS 05506 seen via CO mapping harbors 
several hundreds of solar mass of material (for D=2.8kpc) 
within a region 3 times bigger in size than that probed 
by our ammonia data. Although CO masses are not the 
best measure of the mass of a cloud core, the total mass of 
the associated cloud is likely well above what is needed to 
support massive star formation. It is perhaps unwarranted 
to be too concerned about the relatively low value of the 
cloud core mass, given the possibility that the core has 
been strongly altered or dissipated in the star formation 
process, especially in the presence of the observed strong 
outflows. 

The multiple optical knots in IRAS 05506 imply that 
material has been ejected from the central source in the 
form of high-speed bullets rather than in a continuous 
fashion. The presence of this wide-angle spray of bullets 
in IRAS 05506 indicates, by analogy with Orion, that the 
launching source may be a massive protostar. A massive 
star, depending on its evolutionary status, could produce 
an ultra-compact HII region with detectable continuum 
emission. But since we do not find a measurable radio 
continuum towards IRAS 05506, either it is not a massive 
star, or (i) the central HII region is still very small, (ii) the 
source is sufficiently distant that its radio continuum is be- 
low our noise limit, (iii) the central star has been caught 
at a very early age while its temperature is still too low to 
provide sufficient UV flux for ionisation - this possibility 
is consistent with the late spectral type of the central star 
inferred from our optical spectroscopy. 

We compare IRAS 05506 with the high-mass protostel- 
lar objects (HMPOs) identifled by Sridharan et al. (2002, 
hereafter Setal02). These HMPOs are dense molecular 
cores detected in a CS J=2-l survey of massive star- 
forming regions, which are bright at FIR wavelengths 
(-F60 > 90 Jy) and generally have low (or absent) cen- 
timeter radio fluxes. IRAS 05506 satisfles these selec- 
tion criteria for HMPOs, but shows a signiflcantly smaller 
NH3 line-width (AVArif3=0.96km s"^) and kinetic tem- 
perature (12.5 K) than the median values for HMPOs 
{AYnhs = l-9km s^^ and Tkin = 18.5 K: derived, as for 
IRAS 05506, from single-dish ammonia observations) . In 
fact, Tkin for IRAS 05506 is lower than the smallest value 
for the 40 HMPOs in which Tkin was measured. 

Setal02 have made 2-component dust model flts to the 
IRAS fluxes for these HMPOs and flnd cold dust com- 
ponents with Td 30 — 70 K, and hot components with 
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Td ~ 100 — 250 K, i.e. somewhat higher than the corre- 
sponding ones for IRAS 05506. The bulk of the luminos- 
ity for the HMPOs comes from the cool component, as we 
have found for IRAS 05506. The luminosity of IRAS 05506 
is comparable to those found for many objects in the 
HMPO hst of Setal02. However, the associated IRAS 
fluxes (at 60^m; only upper limits to the lOO/im fluxes are 
available for most HMPOs) and thus the derived masses 
of the cool dust component for HMPOs with comparable 
(near) kinematic distances are generally higher than that 
for IRAS 05506, by factors of ~1.5 or more. 

IRAS 05506 's ammonia line- width and gas kinetic 
temperature are even lower than the median values 
(1.5 km s~^ and 16.9 K) for high-mass starless cores (HM- 
SCs: Sridharan et al. 2005), which are apparently less 
evolved than HMPOs. HMPOs, but not HMSCs, have 
MSX Point Source catalog counterparts, indicative of the 
presence of a luminous, hot, central source. IRAS 05506 
may thus be an extreme example of the HMPO class, lying 
at the bottom end of the ammonia line-width and kinetic 
temperature distributions (see, e.g.. Fig. 3 of Sridharan et 
al. 2005) , and thus at the lower end of the mass range for 
these objects. In a plot of HMPO core mass (derived from 
single-dish mm-continuum data) versus FIR luminosities 
(Fig. 6 of Setal02), there are two sources with core masses 
comparable to that derived from the mm continuum for 
IRAS 05506. Given that the latter is based on interfero- 
metric observations, which resolve out flux on scales more 
extended than about 4", whereas the core masses in Se- 
tal02's Fig. 6 are based on single-dish (beamsize 11") ob- 
servations, IRAS05506's mm-continuum core mass likely 
compares well with even higher core-mass sources in Se- 
tal02's Fig. 6, which has 3 more sources with core masses 
less than about 3OM0. We note that the HMPO list of 
Setal02 may well include stars at the upper end of the 
intermediate-mass range (i.e., 8 Mq); for example, we have 
found at least one example of an HMPO from this list, 
IRAS 20126-1-4104, for which recent observations suggest 
a mass of about 7Mq (Cesaroni et al. 2005). 

The optical detection of IRAS 05506 might seem sur- 
prising, considering the general lack of optical detections of 
HMPOs. It is plausible that the outflows from IRAS 05506 
have cleared away enough circumstellar cloud material 
from the near- vicinity of the central star, allowing a de- 
tectable amount of starlight at optical wavelengths to leak 
out and be scattered into our line-of-sight by dust in the 
walls of the cavity cleared out by the outflow. Note that 
this scenario does not require the outflow cavity to be di- 
rected exactly along the line-of-sight. In addition, we do 
not know the inclination and opening angles of the molec- 
ular outflow, and it is quite possible that these are suffi- 
ciently favorable for the above scenario to be feasible. 

The jet-like feature emanating from the central source is 
directed along an axis roughly orthogonal to the average 
position angle for the most prominent knots. This sug- 
gests that either the collimated outflow which created the 
knots has suffered a drastic change in direction, or Sa-jet is 
the result of an independent high- velocity outflow. If the 
latter hypothesis is correct, it could indicate the presence 
of an unresolved binary star associated with source Sa. 
Taking the lowest value (3150Lo) of the total luminosity 
as derived from the SED modelling (§4.2 1, and dividing it 



equally between the two stars in the binary, we flnd that 
these stars lie on the Bernasconi &; Maeder (1996) pre- 
main-sequence tracks for 9Mq stars. Hence, even if the 
central star is a binary, at least one of the components is a 
massive star (i.e., if D is 2.2 kpc or larger). Note that the 
expansion time-scales derived for the fast molecular out- 
flow and the fast-moving optical knots are very similar, in- 
dicating that these were initiated quasi-simultaneously in 
IRAS 05506's recent history. So if the optical and molecu- 
lar outflows are produced from different stars in a binary 
system, their similar expansion time-scales suggest a com- 
mon triggering mechanism for launching them. 

6. CONCLUSIONS 

We have made a serendipitous discovery of an enigmatic 
outflow source, IRAS 05506-f 2414, as part of a multi- 
wavelength survey of pre-planetary nebulae. Subsequent 
to the HST imaging, we carried out optical spectroscopy, 
millimeter- wave interferometry of the CO J=l-0 line, as 
well as water and methanol maser line and ammonia in- 
version line observations at ~ (22 — 25) GHz. 

(1) The HST optical and near-infrared images show a 
bright compact central source, and a fan- like spray of high- 
velocity (projected radial velocity 350km s~^) elongated 
knots which appear to emanate from it. These struc- 
tures appear to be analogous to the near-IR "bullets" 
around the Orion KL nebula and to have similar masses 

1O~^M0). The central source shows a jet-like exten- 
sion along PA~230°, which is roughly orthogonal to the 
average direction of the bullet spray seen in the optical 
images. 

(2) The high-resolution (3'.'9x2'.'9) CO J=l-0 data show a 
compact emission core of size (< 5"); the line profile has 
a narrow core and weak, high- velocity wings. The latter 
represent emission from a high- velocity bipolar outflow di- 
rected along the optical jet feature. A 10 mJy continuum 
source (at A = 2.6 mm) is also seen at the same location 
as the line peak. The continuum source structure consists 
of an unresolved core and a faint extension whose orien- 
tation is consistent with the optical jet and high-velocity 
molecular outflow direction. 

(3) H2O maser emission was detected both with the VLA 
and the GET. With the latter, NH3(1,1) and (2,2) lines 
were detected, but no methanol emission was found (rms 
antenna temperature was ~25mK). The H2O maser emis- 
sion is spatially unresolved. Our map of NH3 emission 
shows an extended molecular cloud of size about 1', with a 
kinetic temperature of 12.5 K, and a mass of about IIM©, 
assuming a fractional NH3/H2 abundance ratio (by num- 
ber), f{NH3) = 4 X 10-^ and a distance of 2.8 kpc. In 
comparison, the virial mass of the cloud is 48(73) Mq, for 
a r^^ (r^^) density distribution. The factor 4 (or more) 
discrepancy between the virial cloud mass and that de- 
rived from the ammonia data may be attributed to the 
large uncertainty in our assumed value of /{NH^) and the 
distance, and/or to non-gravitational linewidths, such as 
unresolved outflow motion. 

(4) The 8 GHz continuum observations show no contin- 
uum emission (3cr upper limit of 150 /iJy) within a fleld 
of approximately 40" diameter centered on the position of 
IRAS 05506. 

(5) IRAS 05506 is most hkely not a PPN. The combined 
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presence of H2O and OH maser, far-IR, NH3, CO, CS and 
C^^S emission strongly argues for a dense, star-forming 
core being associated with IRAS 05506. IRAS 05506 is 
most likely either a massive or intermediate-mass YSO, 
but not a low-mass YSO. If it is a massive star then 
the lack of radio continuum and the late G - early K 
spectral type we find for the central star implies that in 
IRAS 05506, we are witnessing the earliest stages of a mas- 
sive star's life, while its temperature is still too low to 
provide sufficient UV flux for ionisation. 

(6) IRAS 05506 resembles objects in a recently identified 
class of high-mass protostellar objects (HMPOs). Given 
the relatively low gas kinetic temperature and narrow line- 
width derived from the ammonia data, IRAS05506 could 
be an extreme member of this class, lying at the lower end 
of the mass range for these objects. 

(7) A multi-component dust model provides a reasonable 
fit to the mid-to-far infrared SED of IRAS 05506; we find 
that the bulk of the mass is in a cool dust component 
with a temperature of 27 K and a total (gas-|-dust) mass 
of 26 M0(and about 1% of this in a warmer component at 
79 K), assuming a gas-to-dust ratio of 100. The 2.6 mm 
continuum flux measured with OVRO gives a total mass 
of 2.2-6.4 MqIu a ~4"-size region around source Sa, as- 
suming dust temperatures in the range 79-27 K. 

(8) Both the molecular high- velocity outflow and that pro- 
ducing the optical knots have extremely short and similar 
expansion time-scales of about 200 yr, suggesting that ei- 
ther these are manifestations of a common ejection event 
in IRAS05506's recent past, or if they are due to indepen- 



dent ejections, there must have been a common triggering 
mechanism for launching them. The molecular outflow 
cannot be driven by radiation pressure since the expan- 
sion time-scale is orders of magnitude smaller than the 
time required by radiation pressure to accelerate the ob- 
served bipolar outflow to its current speed. 
(9) We have estimated physical parameters of the high- 
velocity bipolar outflow from the interferometric CO data 
~ our estimates are lower limits because the CO emis- 
sion is optically thick. We find that the outflow mass 
is, Mbip > O.I8M0, the mass-outflow rate is, M^p > 
9.7 X 1O~'*M0 yr~^, the scalar momentum is, Pgc > 
4.3M0km s""'^, the mechanical energy is, Emech > 4.6 x 
j^q45 gj-g^ g^j-^^ |-]^g mechanical power is, L^ech > 7.9 x 10^^ 
erg s~^. 
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Fig. 1. — HST/ACS images of IRAS 05506. Panels a & b show the optical images taken with ACS {log stretch). The position (J2000 RA 
= 05:53:43.54, Dec = 24:14:45.2) of the H2O maser emission feature detected with the VLA is shown as a small circle in panel a; the size of 
the circle denotes the positional uncertainty. An expanded view of source Sf, (inset, panel b) shows a close binary pair with a separation of 
0^'12. Panel c shows the FllOW taJien with NICMOS {log stretch) — a small square patch to the right of Sb corrupted by an instrumental 
artifact has been blanked out. A few of the features in the panel a are also labelled in c for visual registration. Panel d shows expanded views 
(0'.'63x0'.'63) of the knots Kl - KA extracted from the F606W image {linear stretch). 




Fig. 2.— FllOW NICMOS image of IRAS 05506 highlighting the presence of knotty structures (fclO - fcl4) not seen in the F606W ACS 
image (some of the common features {Kl, K2, Sa, and Sb are also labelled for visual registration). The greyscale shows log{log(I)), where / is 
the intensity; above a threshold value, the color stretch cycles repeatedly through the range of colors utilized, giving the mottled appearance to 
the brightest regions in the image. Faint stars in this image are labelled /si — fsA. The dashed circle covers a central region around Sa where 
knotty structures cannot be reliably distinguished from PSF structures. The dotted arcs delineate angular wedges where faint filamentary 
structures can be seen; many of these appear to be radially directed away from Sa. The spiky linear structures oriented at ±45° to the 
horizontal (denoted by dashed black lines) which appear to emanate from Sa are diffraction spikes produced by the telescope optics. Note 
that the orientation of this image is different from the one in Fig.^ 



12 





Fig. 4. — Extended nebulosity in the vicinity of IRAS 05506 is seen in this color-composite image made using 2MASS J (blue), H (green) 
& Ks (red) band images. Each individual image is displayed on a log stretch). Compact sources Sa and Sj, are marked for reference. 




Fig. 5. — (a) FllOW and (b) F160W NICMOS images of IRAS 05506 (shown with a square-root stretch and false-color intensity scale) 
highlighting the presence of low-contrast features. The intensity in the very bright central region surrounding Sa has been scaled down by a 
factor 0.06 (0.03) in the FllOW (F160W) image in order to better show the structure of the surrounding faint nebulosity. The large ellipse 
(panel a) encircles a tube-like structure (Tube); inset in panel b (lower right corner) covers a 2.5" X 4.3" patch around the Tube with a 
different stretch to show its faint parts more prominently (inset shows I^'^^, where / is the intensity). The dashed straight line (panel a) 
marks a roughly linear structure, L, joining two locally bright regions (Bl and B2). The small ellipse in panel a shows the location of B2, 
which is more diffuse and extended compared to Bl. The spiky linear structures oriented at ±45° to the horizontal which appear to emanate 
from Sa, including the 3 bright yellow/orange streaks seen to the east of Bl, are diffraction spikes produced by the telescope optics. Note 
that the orientation of these images is different from the one in Fig.[l] 
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Fig. 6. — Keck/ESI long-slit spectra of IRAS 05506 showing high-velocity outflows and bow-shock structures in the Ha and [Nil] lines. The 
narrow vertical features are due to night-sky emission lines. 
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Fig. 7. — Average intensity spectrum extracted from a 3" x 3" box centered on the peak intensity position in the OVRO CO J=l-0 map. 
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Fig. 8. — OVRO CO J=l-0 (contours) and 2.6 mm continuum (grey-scale) map of IRAS 05506 (5a). The clean beam for the CO (continuum) 
map has a FWHM of 3'.'9 X 2^9, PA -24° ( 4^'6 X 3"2, PA -15°). The yellow, red and cyan contours show the CO intensity integrated over 
the Visr velocity ranges [2.6,10.4], [49.4:119.6], and [-36.4:-106.6] km s~^, respectively. The corresponding values of the peak contour intensity 
and (spacing) are 6.3 (0.66), 3.8 (0.65), and 4.8 (0.65) (Jy/beam)km 8"^ 
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Fig. 9. — Peak H2O emission spectrum of IRAS 05506 observed with the GBT. 
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Fig. 10. — The ammonia emission observed with the GBT at the peak raster position, (a) The (2,2) transition, boxcar smoothed by a factor 
4 (to a resolution of 0.32 km s~^), and resampled by the same factor, (b) The (1,1) transition, without smoothing (resolution of 0.08km s"'^). 
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Fig. 11. — ISO SWS and LWS spectra of IRAS 05506, showing a deep 10 /^m amorphous silicate absorption feature, and emission in the 
OI (at 63 and 146 /im; dashed vertical bars) and CII (158 /im: dashed vertical bar) lines as well as numerous far-infrared CO, OH and H2O 
lines. The locations of the CO lines from J=14— 13 (186.0/im) to J=26-25 (100.5/.im) are shown by short (solid) vertical red bars, and those 
of the most prominent lines of H2O (both ortho- and para-) and OH are marked by green (dotted) and blue (dash-dot-dot) vortical bars, 
resprectively. The symbols show broad-band photometric data from IRAS. The red curve shows a simple dust model fitted to the data. 



